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Abstract 

Recently, Y. Sobouti (2007) has provided a metric theory f(R) that can ac- 
count for certain dynamical anomalies observed in spiral galaxies. Mendoza & 
Rosas-Guevara (2007) have shown that in this theory there is an extra-bending as 
compared to standard general relativity. In the present work we have developed in 
more specific detail this additional lensing effect and we have made evaluations of 
the a parameter used in the model adjusting the t heory to observa tions in X-rays 
of 13 clusters of galaxies with gravitational lensing ([Hoekstral ((2007J)). 



1 Introduction 

Some modified theories of gravity are an actual alternative to the dark matter paradigm, 
which in spite of its success to explain astrophysical problems like the flat rotation curves 
of spiral galaxies, the gravitational lensing by clusters of galaxies and the formation of 
structure in the early Universe, is not the ultimate word of the story. 

In the different alternatives we find the modifications to the Hilbert-Einstein action, 
which introduce an arbitrary function of the Ricci scalar R in the action. This modification 
is relevant, because there's no exist an a priori reason to restrict the gravitational action 
to R and has the advantage to be a metric theory of gravitation that carry out with several 
physical principles and at the same time accounts to the energy-momentum conservation, 
as well as the equations of the theory are invariants. 

In this sense, Sobouti (2007) has developed a metric theory f(R) that reproduces flat 
rotation curves in different spiral galaxies (with f(R) = R 1 ~ a l 2 and a <C 1), reproduces 
the Tully-Fisher relation, converges to a version of MOND in the weak field approximation 
and the resulting differential field equations are of the second order. 
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2 Modified field equations 



he m o del introduces the follow ing; modification to the Hilbert-Einstein action fjCapozziello 



(120021 ): ICapozziello et all (120031 )): 



s 9 = -J~f(R)V=gdn. 



Variation of the action S = S n 



field equations (ICapozziello et al.l (120031 )): 



S m with respect to the metric g^ u gives the following 



f'(R)R»„ - \g, v = T, v + (f(R)).^ - (/' (R)) ;X V- 



(2) 



For an application of this theor y to different g a 



a Schwarzschild-like one given by (jCognola et al.l (120051 ) 



actic s ystems, the me tric is chosen like 



Sobout: 



3j2Qo3)): 



ds 2 = B{r)dt 2 - A(r)dr 2 - r 2 (d6 2 + sm 2 6dip 2 ) . 



(3) 



If the function / differs not too much from general relativity at first order of approx- 
imation, the exterior solution (p = = p) to the field equations is obtained with the 
assumption that v = df/dR is a parametric function v(r,a) with a a real number such 
that a < 1 (for a = the solution converges to G.R.). Proposing a solution of the form 
v oc r@,P G R, results in v ~ r a and the solutions for the metric are (jSoboutil (120071 )): 



Mr) 



1 — a 



r 



Bfr) 



A(r)' 



where s := 2GM is the Schwarzschild radius. For a < 1, f(R) is given by: 



(4) 



f(R) = {3a) a/2 s- a R l - a/2 ~Ry_- ^Ln (s 2 R) + ^Ln(3 



a 



(5) 



3 The deflection of light and gravitational lensing 

Mendoza & Rosas-Guevara(2007) have shown that the deflection angle f3 for a spherically 
symmetric mass distribution in this metric theory (cf. eq.3), is given by: 



P = 7T 



2yT^_ 

2 - a 



a 



+ 2^/1 -a 



l-a/2 



(6) 
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where r m is the impact parameter (the closest approach to the lens). When a = 0, the 
deflection angle equals the Einstein's one: (3e = 2s /r m . 

For an application in gravitational lensing for a symmetric object (e.g. a c l uster of 
galaxies) with mass M, the lens equation is given by (IMendoza fc Rosas- Guevara (120071 )): 



e 2 -e($ + c 1 )-c , 2 e a/2 = o, 



(7) 



where: G\ :- 



Dr 



e E d s 



2-a 



4GM D LS . 
D S D L ' 



6 := 6/6e] $ := i/8e', and 7 is the angular position of the source, 9 is the angular 
position of the image, D$, Dl and Dl$ are the observer-source, observer-lens and lens- 
source distances, respectively. 

The solution to the lens equation at first order, with Oi a small linear perturbation to 



the so lution of general relativity, G>o = \ ($ ± V^ 2 + 4) , is given by (IMendoza &: Rosas- Guevara 
( 120070 ): 



ei 



a/2 



29 - $ - Ci - f CsGn 72 " 1 ' 



(8) 



With this perturbation the magnification factor is (IMendoza fc Rosas-Guevaral (120071 )): 



fi = det [J(0)] _1 



6d6 
$d¥ 



e + de 



(9) 



$ d$ ' 

where J (9) is the jacobian of the transformation 9 — > 7 (from the lens plane to the source 
plane) . 

For det [J(0)] = (when the magnification factor /1 diverges), the critical curve for 
$ = and for 0q~(O) = 9 E , is given by: 
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and for 0jT (0) = 2vr - 9 E by: 



c 



2 (e -) 2 -i + c 2 (e )^ (i-f) 



1 + c 2 (1 - f ) 

2 - c x - f c 2 



0, 



(10) 



-\a/2 



9. 



c 



9i 



[ID 



2e --c 1 -fc 2 (eo) a/2 1 

From the symmetry of the problem we have (27r — 9 C ) = 9^ := 0c, and we can take 
9c as the equivalent to the Einstein angle 9e, i-e. an "Einstein ring" in this theory have 
an aperture of 9c, and 9c = 9e for a = 0. 
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4 Approximation method to a 



Since these results are derived from a metric with spherical symmetry, we use observations 
in X-rays of gr avitationa l lensin g by 13 clusters of galaxies that approximate well to this 



symmetry (see 



Hoekstra 



(120071 ) and references therein). In order to estimate the param- 
eter a, we assume that the estimated values of the Einstein angle from the observations 
approximates 9c- We use this as a first way to mesure the order of magnitude to a. 

The cosmological model used is: H = 70kms~ 1 Mpc~ 1 , Q m = 0.3 and f^A = 0.7. To 
calculate the Einstein angle we use the virial mass of the sy s tem w i th the mean velocity 



disper sion observed in a particular system (IBorgani et al.l (119991 ); 
( I200lk 



Girardi Sz Mezzetti 



Using this method we can solve numerically the equation (10) to obtain a. 



5 Results 

Sobouti's f(R) theory is a first approximation to a modified metric theory of gravity 
in order to account for phenomena usually adscribed to dark matter. However, there 
is a caveat with Sobouti's approach. The mass dependency of a destroys one of the 
most important facts of an ordinary action: it should be invariant under the changes of 
sources. Also, this mass dependency means that the function f(R) varies from source to 
source. However, as mentioned by Sobouti (2007), perhaps we need to re-think whether 
this postulate has to be accepted when modifications to standard general relativity are 
done. 

The way Sobouti (2007) calculated the value for the parameter a was using a spher- 
ically symmetric metric, applied to spiral galaxies. This is of course wrong since these 
galaxies are far from spherical symmetry. This is the reason as to why we decided to 
callibrate the a parameter using spherical symmetric astrophysical systems: cluster of 
galaxies. From the many cluster of galaxies known we have chosen a set of very well 
spherically symmetric objects. 

To summarize, our main results are as follows: 

• The caustic in the source plane is the point colinear with the lens and the observer. 
This is an expected result since the metric is a Schwarzschild-like one for a spherically 
symmetric matter distribution. 
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Figure 1: Plot of a vs. (M/1O 14 M ) 1 / 2 . The best linear fit is shown by the line. 

• The a-mass relation from Sobouti (2007) is given by a = a ( J 

Our best linear fit gives a value of «o — (3-5 ± 0.3)xl0~ 9 (see Fig. 1), that differs 
in three orders of magnitude from the value obtained by Sobouti (2007), cto = 
(3.07 ± 0.18)xl0~ 12 , for 31 spiral galaxies with a = Iv^jc? {y^ is the asymptotic 
constant tangential speed of the galaxy) . 

• The log-log plot of a vs. M/M , leads the power law fit: a = [(1.47 ± 0.51)xl0~ 5 ] 
i.e. a relation of the form: 



0.25 



a(M) = ai 



M 



1/4 



(12) 



Although the calculations are only a first approximation to the order of magnitude 
to the a parameter, this is the correct way for evaluating it, since the distribution of 
barionic matter approximates to a spherical symmetric one in clusters of galaxies. 

One could be tempted to use the results also for quite "spherical" elliptical galaxies. 
However, this approximation fails because the Einstein angle in these galaxies is 
very small and so the variables C\ and Ci (cf. eq.7) diverge. 
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Figure 2: Plot of log(a) vs. log(M/M ). The best linear fit is shown by the line. 

• Future work : In section H] we approximated the value a by using the estimated 
Einstein angles and equate them to Qq- This is a first approximation to the real 
value of a and the only way to correct this is to develope a complete theory of 
gravitational lensing for the metric theory proposed by Sobouti. In this way, we 
intend to determine with much precision the exact value for the a parameter, and 
check whether this theory can account for different astrophysical problems. 
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